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£h Abstract 

h— 5 Given a collection of n opaque unit disks in the plane, we want to stack them on top of 

each other in an order so as to maximize the total length of all pieces of their boundaries 

visible from above. We call this total length the visible perimeter of the collection with 
respect to the particular order. We prove that if the centers of the disks form a dense 
point set, i.e., the ratio of the maximum distance to the minimum distance between them 
is (^(n 1 / 2 ), then there is a stacking order for which the visible perimeter is f2(n 2 / 3 ). We 
c/3 also show that the order of magnitude of this bound cannot be improved in the case of 

the n 1 / 2 by n 1 / 2 piece of a sufficiently small square grid. On the other hand, if the set 
of centers is dense and the maximum distance between them is small, then the visible 

t— I perimeter is 0(n 3 / 4 ) with respect to any stacking order. We construct collections of disks 

for which the order of magnitude of this bound is tight. The above results partially answer 
some questions raised by Cabello, Haverkort, van Kreveld, and Speckmann. 

1 Introduction 

O 

^vq In cartography and data visualization, one often has to place similar copies a symbol, typically 

an opaque disk, on a map or a figure at given locations [De99] , |Gr90] . The size of the symbol is 
J> sometimes proportional to the quantitative data associated with the location. On a cluttered 

map, it is difficult to identify the symbols. Therefore, it has been investigated in several 
studies how to minimize the amount of overlap |GrC78| . [SIMPS] . 

In the present note, we follow the approach of Cabello, Haverkort, van Kreveld, and 
Speckmann [CaHlOj. We assume that the symbols used are opaque circular disks of the same 
size. Given a collection D of n distinct unit disks in the (x, y)-plane, a stacking order is a one- 
to-one assignment / : D — )■ {1, 2, . . . , n}. We consider the integer f(D) to be the z-coordinate 
of the disk D £ T>. The map corresponding to this stacking order is the 2-dimensional view of 
this arrangement from the point at negative infinity of the z-axis (for notational convenience, 
we look at the arrangement from below rather than from above.) In particular, for the lowest 
disk D, we have f{D) = 1, and this disk, including its full perimeter, is visible from below. The 
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Figure 1: Left: A collection of unit disks in the plane. Right: A stacking order for them. 



total length of the boundary pieces of the disks visible from below is the visible perimeter of T> 
with respect to the stacking order /, denoted by visible(P, /). We are interested in finding a 
stacking order for which the visible perimeter of T> is as large as possible. See Figure [l} Given 
an integer n, we define 

v(n) = inf max visible (T>, /), 

\V\=n f 

where the maximum is taken over all stacking orders /. We would like to describe the asymp- 
totic behavior of v(n), as n tends to infinity. 

Cabello et al. have already noted that v(n) = f}(n 1//2 ); in other words, every set T> of 
ra disks of unit radii admits a stacking order with respect to which its visible perimeter is 
Q(n 1//2 ). Indeed, by a well-known result or Erdos and Szekeres [ErSz_35j, we can select a 
sequence of [~ra 1//2 ] disks Di £ T> (1 < i < |~n 1//2 ~|) such that their centers form a monotone 
sequence. More precisely, letting Xj and yi denote the coordinates of the center of Z)j, we have 
x\ < X2 < X3 < . . . and either y\ < y2 < 2/3 < . . . or y\ > 1/2 > 2/3 > . . .. Then, in any stacking 
order / such that f(Di) = i for every i, 1 < i < ["ra 1 / 2 ] , a full quarter of the perimeter of each 
Di (1 < i < [n 1 / 2 ] ) is visible from below. Therefore, the visible perimeter of T> with respect 
to / satisfies 

visible(P,/) > ^\n 1/2 ]. 

At the problem session of EuroCG'll (Morschach, Switzerland), Cabello, Haverkort, van 
Kreveld, and Speckmann asked whether v(n) = 0(n); in other words, does there exist a 
positive constant c such that every set of n unit disks in the plane admits a stacking order, 
with respect to which its visible perimeter is at least era? We answer this question in the 
negative; cf. Theorems 2 and 5 below. 

Given a set of points P in the plane, let T>(P) denote the collection of disks of radius 
1 centered at the elements of P. For any positive real e, let eP stand for a similar copy 
of P, scaled by a factor of e. For a stacking order / of T>(P) we will study the quantity 
visible (2? (eP), /). (Note the slight abuse of notation: We denote the stacking order of T>(P) 
and the corresponding stacking order of T>(eP) by the same symbol /. The two orders are also 
identified in Lemmas 1 and 2.1 and in Theorems 2, 3, and 5.) It is not hard to verify that, as e 
gets smaller, the function visible (P(£-P), /) decreases. To see this, it is enough to observe, as 
was also done by Cabello et al. (unpublished), that as we contract the set of centers, the part 
of the boundary of each unit disk visible from below shrinks. As we will see in Lemma 2.1, 
the limit in the following lemma has a simple alternative geometric interpretation. 
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Lemma 1. For every point set P in the plane, and for every stacking order f of the collection 
of disks V{eP), we have 

visible (DfeP),/) > Iimvisible(D(eP),/). 

£->0 



As in |AlKP89j . [Va92| . and |Va96j . we consider C -dense n-element point sets P, i.e., point 
sets in which the ratio of the maximum distance between two points to the minimum distance 
satisfies 

max(|pg| :p,q£P) < 1/2 
min(|pg| : p, q G P,p ^ q) ~ 

Theorem 2. For any C -dense n-element point set P in the plane and for any stacking order 
f , we have 

limvisible(P(eP),/) < C"n 3/4 , 

£-►0 

where C is a constant depending only on C . 

The order of magnitude of the upper bound in Theorem 2 cannot be improved: 

Theorem 3. For every positive integer n, there exists a 2-dense n-element point set P n in 
the plane and a stacking order f such that 

limvisible(P(eP n ),/) > n 3/4 /2. 



In the general case, where P is an arbitrary n-element point set in the plane, we have been 
unable to improve on the easy lower bound 

maxvisible(£>(P),/) = Q(n 1/2 ), 

sketched above. However, under special assumptions on P, we can do better. 

Theorem 4. Every C -dense n-element point set P in the plane admits a stacking order f 
with 

visible(P(P),/) > C"n 2/3 , 

where C" > depends only on C . 

In particular, Theorem 4 provides an fi(n 2 / 3 ) lower bound for the visible perimeter of a 
collection of n unit disks centered at the points of an n 1 / 2 x n 1 / 2 piece of a square lattice of 
any side length, under a suitable stacking order. Assuming the side length is very small, this 
is better than the line-by-line "lexicographic" stacking order for which the visible perimeter 
is only 0(n 1//2 logn). It turns out that in this case there is no stacking order for which the 
order of the magnitude of the visible perimeter would exceed re 2 / 3 . 

Theorem 5. Let re be a perfect square and let G n denote an n 1 / 2 by n 1 / 2 piece of the square 
lattice in the plane. For any stacking order f , we have 



limvisible(£>(eP n ),/) = 0(n 

£->0 



2/3^ 
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Consequently, we have v(n) = 0(n 2 / 3 ). 

In this note, we are concerned only with worst case results. It turns out that in the worst 
case, the centers of the unit disks are very close to each other, so all disks have a point in 
common. This is, of course, not a realistic assumption in the labeling problem in cartography 
that has motivated our investigations. In practical applications, only a bounded number of 
unit disks share a point. In the last section, we show how to apply our results when such a 
bound is known. 

In Section 2 we establish Theorems 2 and 3. The proof of Theorem 4 is presented in 
Section 3. In Section 4 we consider the square grid and present a much simpler proof of this 
special case of Theorem 4 based on Jarnik's theorem |Ja25j : we then prove Theorem 5, which 
states that the bound of Theorem 4 is tight in this case. The last section contains concluding 
remarks and open problems. 

2 Dense sets with largest visible perimeter 
Proofs of Theorems 2 and 3 

First, we express the limit of visible perimeters appearing in Theorems 2 and 3 in a simpler 
form. Given a set of points P in the plane, let convP stand for its convex hull. Let D(p) 
denote the unit disk centered at p and let T>(P) stand for the set {D(p) : p G P}. 

Fix an orthogonal system of coordinates in the plane. For any point p = (x, y) and for any 
e > 0, let ep denote the point with coordinates (ex,ey). 

Lemma 2.1. Let P = {p\,p2, ■ ■ ■ ,p n } be a set of points in the plane, let e > 0, and let f 
be the stacking order of the unit disks about the elements of eP such that f(D(epi)) = i for 
i = 1, 2, . . . , n. 
We have 

n 

limvisible(P(e J P),/) = V n, 

E— >0 — ' 

8=1 

where t\ = 2ir, and for all other indices, t% = if pi belongs to convjpi , £>2 > • • • >Pi-i}, and T{ 
is equal to the external angle of the convex polygon conv{pi,p2, • • • ,Pi} at vertex pi, otherwise. 
See Figure^ 

Proof. We prove that the contribution of T>(epi) to the visible perimeter tends to Tj as e — > 
for each 1 < i < n. 

Since D{ep\) is the lowest disk in T>(eP), its whole boundary is visible from below. There- 
fore, its contribution is 2ir. Let i > 1. If belongs to the interior of conv{pi,p2 ; • • • ,Pi-i}, 
then there is a threshold Eo > such that 

i-l 

D(epi) C |J D(e Pj ), 
j=i 

for every e < £q. In this case, no portion of the boundary of D(epi) is visible from below, 
provided that e is sufficiently small. If pi lies on the boundary of convjpi , p2, ■ ■ ■ , Pi}, then it is 
in between some points pj and p^ with 1 < j < k < i and although T>(epi) will not be entirely 
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Pi 




Figure 2: If pi lies outside the convex hull of the preceding points, then n is defined as the 
external angle of the polygon convjpi, . . . ,pi} at vertex p^. 




Figure 3: The triangle Piqq' lies entirely outside the convex hull of p±, . . . ,Pi-i- 

covered by earlier disks for any e > 0, the part of its boundary outside V(epj) UV(epk) tends 
to zero as e — > 0. 

Finally, if pi lies outside convjpi, . . . ,Pi-i}, then it is a vertex of convjpi, . . . ,pi}. Consider 
the external unit normal vectors to the two sides of convjpi, . . .p{\ that meet at epi (or in 
the case if the convex hull is a single segment, the two unit normal vectors for this segment). 
Drawing these vectors from epi, the arc on the boundary of V{pi) between them is of length Tj 
and it is not covered by U}=i D{spj). Thus, it is visible from below, and, as e — > 0, the total 
contribution of the remaining part of the boundary of D(epi) to the visible perimeter tends 
to 0. □ 

Proof of Theorem 2. Consider a C-dense point set P in the plane and let / be a stacking 
order for V{P). Suppose without loss of generality that the minimum distance between the 
elements of P is 1 and the maximum distance (diameter) is at most Cn 1 / 2 . We write P = 
{piiP2-, ■ ■ ■ ,Pn} with f(V(pi)) = i. Using Lemma 2.1, it is enough to prove Y17=i r * — C'v?^ 
for the angles Tj defined in the lemma. As Tj = whenever pi is contained in convjpi, . . . ,Pi-i}, 
we can assume this is never the case. 

For every i (1 < i < n), let per(i) denote the perimeter of conv{pi,p2, ■ ■ ■ ,Pi}- We define 
the perimeter of a segment to be twice its length and the perimeter of a point to be 0. Let 
2 < i < n, consider the two sides of the polygon conv{pi,p2, ■ ■ ■ ,Pi} meeting at pi, and denote 
by q and q' the points on these sides at unit distance from pi. Since no point of P is closer 
to pi than 1, the triangle Piqq' does not contain any element of {pi,P2, ■ • • ,Pi-i}- Hence, 
conv{pi,p2> • • • ,Pi~i} is contained in the convex region obtained from conv{pi,p2, ■ ■ ■ ,Pi} by 
cutting off the triangle piqq' . (In the degenerate case when convjpi, . . . ,pi} is a segment, 
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q = q' is the point of this segment at unit distance from pi , and the empty "triangle" becomes 
just a unit segment.) This observation implies that the perimeter of conv{pi,£>2 5 ■ ■ • ,Pi-i} 
satisfies 

n t 2 

per(i - 1) < per(i) - \piq\ - \piq'\ + \qq'\ = per(z) - 2 + 2 cos ^ < per(i) 

Here we used that the external angle of the triangle piqq' at vertex p, L is Tj. 
Thus, we have 

r 2 

per(i) — per(i — 1) > 

5 

for all i > 1. Adding up these inequalities, we obtain 

n 

per(n) > 

Since per(n) is at most tt times the diameter of P, that is, per(n) < irCn 1 ^ 2 , we have 

i=2 

Applying the relationship between the arithmetic and quadratic means, we can conclude that 

n / n \ V 2 

i=2 \i=2 / 

Taking into account that T\ = 2ir, the theorem follows by Lemma 2.1. □ 

Proof of Theorem 3. Suppose for simplicity that n = 2k 2 for a sufficiently large positive 
integer k. Let pi,P2, ■ ■ ■ ,P2k be the cyclic enumeration of the vertices of a regular 2fc-gon of 
unit side length centered at the origin O, and let C\ denote its circumscribing circle. Take k 
homothetic copies of this 2fc-gon such that their vertices lie on concentric circles C\, C%, . . . , Ck, 
where the radius of C a+ \ satisfies r{C a +\) = r(C a ) + 1 for every a (1 < a < k). The vertices of 
the 2fc-gon inscribed in C a +\ are denoted by P2ka+i,P2ka+2, ■ ■ ■ ,P2ka+2k, in cyclic order. Let 
Pn = {Pi '■ 1 < i < n }- Fix the stacking order / of V{P n ) so that f(D(pi)) = i for every 
i (1 < i < n). 

The minimum distance between the points of P n is 1, and the maximum distance is at 
most 

2r(C k ) = 2(r(d) + k - 1) = — — - + 2k - 2 < 2^2k = 2n 1 / 2 . 

sm(7r/2K) 

Thus, P n is a 2- dense set. 

From any point p on Ck, draw a tangent line to Ck-i and denote the point of tangency by 
p' . Let r be defined as the angle pOp', so that we have 



r(C fc -i) ^ n 1 / 2 - 1 _ 1/4 

r = arccos — — . > arccos -7- — > n 1 . 

r(C k ) n 1 ' 2 
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Figure 4: Left: Partition of D into annular sectors centered at a point p £ K. Right: Point 
Pi contributes external angle at least j3. 

It remains to estimate the angles Tj that appear in Lemma 2.1. Half of the points pi are of 
the form p. L = P2ak+b with < a < k and 1 < b < k. Consider such a point pi S C a +i with 
a > 0, and draw through it a tangent line to C a . Denote the point of tangency by p\. Notice 
that Tj is at least as large as the angle PiOp'^ which is larger than or equal to r. For a = 0, 
we have t% > tt/2 > r. Thus, for all n/2 points pi of the above form, we have Tj > r > n~ 1//4 , 
which yields Y27=i r « — re3//4 /2- Lemma 2.1 finishes the proof. □ 

3 All dense sets have good stacking orders 
Proof of Theorem 4 

Throughout this section, let P be a C-dense n-point set in the plane. We will define a 
stacking order / for T>(P) for which the external angles Tj defined in Lemma 2.1 satisfy 
Ya=1 Ti — C"n 2 / 3 , for some constant C" > depending only on C. Then the theorem follow 
from Lemma 2.1. 

Assume without loss of generality that the minimum distance in P is 1. Then, since P is 
C-dense, there exists a disk of radius Cn 1 / 2 that contains all of P. Let D be such a disk, and 
let K be a circle of radius 2Cn 1 / 2 concentric with D. 

Given a point p G K, we define a family F = F(p) of annular sectors that disjointly cover 
the plane, as follows: For each positive integer i, let Ki = Ki{p) be a circle centered at p with 
radius in~ 1//6 ; then divide each annulus between two consecutive circles into sectors of angular 
length a = C*n -1//3 for a large enough constant C* (as will be specified below). See Figure [4] 
(left). 

Note that each annular sector that intersects D has area 0(1) (since its radius is ©(n 1 / 2 )). 
Since D has area 0(n), there must be 0(n) sectors that intersect D. Call a sector occupied if 
it contains at least one point of P. 

Lemma 3.1. There exists a point p £ K for which £l(n) sectors of F(p) are occupied. 

Proof. Choose p uniformly at random on K and construct the sectors using p and dividing 
the annuli into the correct length sectors in an arbitrary way. For each point pi G P, define 
the random variable n{pi) to be the number of points of P contained in the sector of F{p) 
that contains pi. We claim that the expected value E\n{pi)\ of n(pi) satisfies 

E[n( Pi )} < k 
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for some constant k. 

Indeed, let R = R Pi (p) be a rectangle centered at pi, with dimensions (/c'n 1 / 6 ) x (A;'n~ 1//6 ), 
and with short sides parallel to the line ppi, for an appropriate constant k'. If k' is large 
enough (but constant with respect to re), then R completely contains the sector of F(p) that 
contains p, t . Thus, it suffices to bound the expected number of points of P in R. Note that, 
as p rotates around K, R rotates around its center together with p. 

Partition the plane into annuli centered at pi by tracing circles around pi of radii 1, 2, 4, 8, . . .. 
The annulus with inner radius r and outer radius 2r contains at most k%r 2 points of P, for 
some constant k2- Each such point has probability at most k^n^ 1 ^^ 1 of falling in R (over 
the choice of p), for another constant k^; therefore, the expected contribution of this annulus 
to n(pi) is at most A^A^rn -1 / 6 . Summing up for all annuli with inner radius r < A/n 1 / 6 , we 
obtain that E[n(pi)] < k for some constant k, as claimed. 

Now, call point pi isolated if n{pi) < 2k. By Markov's inequality, each point pi has 
probability at least 1/2 of being isolated. Therefore, the expected number of isolated points 
is at least n/2. There must exist a p that achieves this expectation, and for it we obtain at 
least n/(4fc) occupied sectors. □ 

Proof of Theorem 4. Fix a point p for which F(p) has Cl(n) occupied sectors. Color the 
sectors with four colors, using colors 1 and 2 alternatingly on the odd- numbered annuli and 
colors 3 and 4 alternatingly on the even- numbered annuli. 

There must be a color for which Q(n) sectors are occupied. Consider only the occupied 
sectors with this color. Let these sectors be Si, S2, . . . , S m , listed by increasing distance from 
p, and for each fixed distance, in clockwise order around p. Select one point pi E P n Si from 
each of these sectors. Let the stacking order / start with these points, that is, f(T>(pi)) = i 
for i = 1, . . . , m. The order of the remaining points in P is arbitrary. 

We claim that each selected point pi contributes an external angle of = f2(re -1 / 3 ), which 
implies that ^ r% = fi(n 2 / 3 ), as desired. 

Indeed, consider the i-th selected point p%. Suppose without loss of generality that p lies 
directly below pi. Let and K^+i be the inner and outer circles bounding the annulus that 
contains pi. Trace rays Z\ and z-i from pi tangent to i£fc_i, touching K^_\ at points q\ and 
(72 • See Figure |4] (right). 

Every point pj, j < i, that is not contained in the same annulus as pi lies below these 
rays. Moreover, the angle j3 that these rays make with the horizontal is 0(re -1 / 3 ): Consider, 
for example, the ray z\. The triangle ppiqi is right-angled, with angle Lpipq\ = j3. We 
have pq\ = ©(n 1 / 2 ) and ppi = pq\ + 0(n -1 / 6 ). It follows that piq\ = 0(n 1 / 6 ), and so 
/3 w tan ft = Piqi/pqi = 0(re -1 / 3 ). 

Now suppose that pi-i lies in the same annulus as pi. If the constant C* in the definition 
of a is chosen large enough, then Pi—i must have a smaller y-coordinate than p^. (In the worst 
case, pi lies near the bottom-left corner of its sector and Pi-± lies near the top-right corner of 
its sector.) 

Thus, pi contributes external angle > (3 = ^(re^ 1 / 3 ), as claimed. □ 
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Figure 5: The triangle PiQi^ is the largest isosceles triangle at point pi that does not intersect 
the interior of Cj_i. 

4 The "worst" dense set: the grid 
Proof of Theorem 5 

In this section, we assume that n is a square number and G n denotes an ra 1 / 2 by n 1 / 2 piece of 
the integer grid. Note that G n is a \/2-dense set consisting of n points. 

As we mentioned in the Introduction, in the special case where P = eG n , Theorem 4 has a 
very simple proof. For T>(eG n ), one can produce a stacking order with large visible perimeter 
using the following greedy algorithm (which can also be applied to any other point set P): 
Set P n = G n , and select a vertex of conv(P n ) whose external angle is maximum. Let this 
vertex be p n , the last element in the desired order f greedy Repeat the same step for the set 
P n -i = P n \ {p n }, and continue in this fashion until the first element p\ gets defined. 

By Jarnik's theorem [Ja25j, every convex polygon has 0(n 1//3 ) vertices in G n . Therefore, at 
each step, the greedy algorithm selects a point pi that makes an external angle r« = ^(n -1 / 3 ). 
Hence, Yl T i = ^(^ 2//3 ) for the order f gvee dy Lemma 2.1 completes the proof. 

Now we turn to the proof of Theorem 5. Our proof is an improved version of the proof of 
Theorem 2. There we were concerned with how the perimeter of the convex hull grows as we 
add the points of our set one by one as prescribed by the stacking order. As is well known, 
the perimeter of a convex set in the plane is the integral of its width in all directions (this 
is known as Cauchy's theorem; see e.g. 1 1 'a A {)■">' , Theorem 16.15). The proof of Theorem 5 
is very similar, but we deal with the the widths in different directions in a non-uniform way. 
The width in a direction close to the direction of a short grid vector is more important in the 
analysis than widths in other directions. 

Proof of Theorem 5. Let G n = {pi, . . . ,p n } be an enumeration of the points of G n according 
to a given stacking order, and let Tj denote the corresponding external angles, as defined in 
Lemma 2.1. According to the lemma, we need to prove that X^=i r * = 0(n 2 / 3 ). Let us 
partition this sum into several parts, and bound the contribution of each part separately. 
Let [n] = {1, . . . , n}. We start with the small angles. Let 

I = {i G [n] | n < n~ 1/3 }. 

Clearly, we have 

^r l <n-n- 1 / 3 = n 2 / 3 . 
iel 

As in the proof of Theorem 2, let C% = convjpi, . . . ,pi} and denote the perimeter of C, by 
per(i). Since G n is a piece of the integer grid, we have per(n) = 4(n 1//2 — 1). Consider only 
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those indices i > 1 that do not belong to Iq. For these indices, we have r< > 0, so that pi 
must be a vertex of Cj. For each such point pj, let Cj denote the smallest number satisfying 
the following condition: the segment connecting the points qi and q[ that lie on the boundary 
of Ci at distance q from p,, intersects Cj-i. (In the case where Cj is a segment, we have 
q% = q'i £ Cj-i-) Note that the segment qiq\ contains a point pj with 1 < j < i. See Figure [5j 
In the proof of Theorem 2, we argued that per(z) — per(z — 1) > r?/5. Now the same 
argument gives that per(i) — per(z — 1) > CjT?/5. Let 

j 1 = {i G H\(J u{l}) \ar i >n~ 1 ^}. 

For i € Ji, we have per(i) — per(i — 1) > Tjre _1 / 6 /5. Using that per(i) is monotone in i, we 
conclude that 

Yj Ti - 5n 1/6 (per(n) - per(l)) < 20n 2/3 . 

ieh 

Let 

/ 2 = [n]\({l}U/ U/i). 

To bound the angles t% for indices i G I2, we need a charging scheme and we need to consider 
the growth of the width of Cj in some specific directions. The width of a planar set in a given 
direction is the diameter of the orthogonal projection of the set to a line in this direction. 
Let us associate the directions in the plane with the points of the unit circle K. We identify 
opposite points of this circle as the widths of the same set in opposite directions are the same. 
This makes the total length of K become n. We define a set of arcs along K as follows. For 
any non-zero grid vector v from the integer grid and for any integer £ > 0, let V v ( denote the 
arc of length symmetric around the direction of v. For any direction a G K, let Pi(a) 
denote the width of Cj in the direction orthogonal to a (i.e., the corresponding projection is 
parallel to a). 

The perimeter per(i) is equal to the integral of Piioi) along the circle K (note that after 
the identification of opposite points the length of K became 7r). We have Pi(a) = pj_i(a), 
unless the direction a is tangent to Cj at the vertex pj. Let f7j denote the arc of directions 
where such a tangency occurs. Clearly, the length of Ui is Tj, and for any arc V that contains 
Ui, we have 

(Pi(a) - pi-i(a))da = per(i) - per(z - 1) > qT?/5. 

For each index i G I2, choose a grid point pj on the segment (Recall that the points 
qi and ^ are at distance Cj from pj, and that there is always a grid point between them.) We 
charge the index i to the pair (f(i),£(i)), where v(i) is the grid vector pointing from pj to pi 
and £(i) is the largest integer such that V v u\m^ contains U{. Notice that \v(i)\ < Ci. 

Note that Ui is symmetric around the direction of the segment qiq\. For the angle /3 
between this direction and the direction of v(i) we have \vi\ sin/? = Cjsin(rj/2) (refer again to 
Figure]^. This implies (3 < Cin/\v(i)\, and hence 2~ 1 ^ < 4/3 + 2ti < 6cfr»/|u(i)|. Finally, we 
also have 

/ (Pi(a) -Pi-i{a))da > arf/5 > 2- < W|u(i)|r i /30. 

Let s(v,£) = where the summation extends over all i G I2 charged to (u, £). The 

integral Jy pi(a)da is monotone in i and grows by at least 2 _ ^|u|rj/30 at every i that is 
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charged to (v,£). We have p\(a) = and p n (a) < (2n) 1//2 , so that the total growth is less 
than 2- £ (2n) 1 / 2 . This implies that s(v,£) < 43n 1 / 2 /|^|. 

Consider the set of all pairs (v, £) such that there is an index i G I2 charged to them. We 
have CjTj < n -1 / 6 , Tj > n -1 / 3 and |v| < q, which implies that |t>| < n 1 / 6 . We proved that 
2~ e < 6ciTi/\v\ < 6n _1 / 6 /|u|. On the other hand, we also have 2~ e > 2n > 2n -1 / 3 . Thus, 
for any given grid vector v, there are at most log(6n 1 / 6 /|f |) possible values of £, where log 
denotes the binary logarithm. 

Hence, 

^ Ti = ^ s (V)< log(6n 1/6 /kl)43n 1 / 2 /bl- 

ieh v,i \v\<n 1 / 6 

Evaluating this sum on the grid vectors v of the integer grid, we obtain 

5> = 0(n 2 / 3 ). 

We have 

£ Tt = T1 + « + E * + E ^ = 27r + °(^ 2/3 ) + °(" 2/3 ) + °(^ 2/3 ) = °( n2/3 )' 

1=1 ie/o iei"i 

completing the proof of the theorem. □ 

5 Concluding remarks 

A. The greedy algorithm described at the beginning of Section 4 was first considered by Cabello 
et al. (unpublished) in the context of maximizing the minimum visible perimeter of a single 
disk. They showed that the order /greedy is always optimal for this purpose. Unfortunately, 
this stacking order is not always optimal with respect to the total visible perimeter. Indeed, 
let n be a perfect square and consider the set of points {pi \ 1 < i < n}, where the polar 
coordinates of pi are (e bl , 27ri/n 1 / 2 ) with b > sufficiently small. This point set is obtained 
as the intersection of n 1 / 2 equally-spaced rays emanating from the origin, and n 1//2 "rounds" 
of a very tight logarithmic spiral centered at the origin. The greedy algorithm produces the 
stacking order indicated by the indices, so it takes the points of P outwards along the spiral. 
The contribution T{ is equal for every point pi with i > n 1 / 2 and tends to 27r/n 1//2 as b goes 
to zero, making X] r « = ©( n1 ^ 2 ) if b is small enough. However, taking the points ray by ray 
in a cyclic order, going outwards along each ray, the contribution T{ is a constant for the first 
half of the points, making J2 T i = 

B. Theorem 4 can be generalized to point sets satisfying weaker density conditions. Indeed, let 
P be a set of n points in the plane with diameter D and minimum distance d. A randomized 
construction, similar to the one used in the proof Theorem 4, guarantees the existence of a 
stacking order / such that visible(X>(P), /) = Q(n/(D/d) 2/3 ). This beats the ftfn 1 / 2 ) bound 
mentioned in the Introduction as long as D/d = o(n 3 / 4 ). 

C. We have seen that in all collections V of n disks, for which max j visible (V, f) is close 
to being minimal (close to v(n)), all disks share a point. Of course, this is not a realistic 
setting for applications. Nevertheless, studying the behavior of the function v(n) (defined in 
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the Introduction) also helps in minimizing max f visible (T>, f) under the assumption that no 
more than c disks in T> have a point in common. In this case, there exists a stacking order / 
with 

visible(£>, /) = n(v(c)n/c). 

To see this, place a random translate of an infinite square grid of side-length 4 on the 
plane. For each unit disk, the probability that it belongs entirely to a single cell is at least a 
positive constant. Thus, we can fix the grid in such a way that the number of disks that lie 
entirely in one cell is f2(n). If a cell contains k disks, we have k = 0(c), and the best stacking 
order makes their visible perimeter at least v(k) = Q(kv(c)/c). Starting with an optimal order 
of these disks and following with an arbitrary order of the remaining ones, the claimed bound 
follows. 

This bound is tight, as shown by the union of n/c "worst case" sets of c disks far from 
each other. 

Acknowledgement. The authors express their gratitude to Radoslav Fulek and Andres Ruiz 
Vargas (EPFL), for many useful discussions on the subject. 
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